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ABSTRACT 

Polaris,  the  nearest  and  brightest  classical  Cepheid,  is  a  single-lined  spectroscopic  binary  with  an  orbital  period  of 
30  yr.  Using  the  High  Resolution  Channel  of  the  Advanced  Camera  for  Surveys  on  board  the  Hubble  Space  Telescope 
( HST)  at  a  wavelength  of  ~2255  A,  we  have  directly  detected  the  faint  companion  at  a  separation  of  0"17.  A  second 
HST  observation  1 .04  yr  later  confirms  orbital  motion  in  a  retrograde  direction.  By  combining  our  two  measures  with 
the  spectroscopic  orbit  of  Kamper  and  an  analysis  of  the  Hipparcos  and  FK5  proper  motions  by  Wielen  et  al.,  we 
find  a  mass  for  Polaris  Aa  of  4.5+2, \  M© — the  first  purely  dynamical  mass  determined  for  any  Cepheid.  For  the  faint 
companion  Polaris  Ab  we  find  a  dynamical  mass  of  1.26+ q^4  Mq,  consistent  with  an  inferred  spectral  type  of  F6  V 
and  with  a  flux  difference  of  5.4  mag  observed  at  2255  A.  The  magnitude  difference  at  the  V  band  is  estimated  to  be 
7.2  mag.  Continued  HST  observations  will  significantly  reduce  the  mass  errors,  which  are  currently  still  too  large  to 
provide  critical  constraints  on  the  roles  of  convective  overshoot,  mass  loss,  rotation,  and  opacities  in  the  evolution 
of  intermediate-mass  stars.  Our  astrometry,  combined  with  two  centuries  of  archival  measurements,  also  confirms 
that  the  well-known,  more  distant  (18")  visual  companion,  Polaris  B,  has  a  nearly  common  proper  motion  with  that 
of  the  Aa,  Ab  pair.  This  is  consistent  with  orbital  motion  in  a  long-period  bound  system.  The  ultraviolet  brightness 
of  Polaris  B  is  in  accordance  with  its  known  F3  V  spectral  type  if  it  has  the  same  distance  as  Polaris  Aa,  Ab. 
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1.  INTRODUCTION 

Cepheid  variable  stars  are  of  central  importance  in  galactic 
and  extragalactic  astronomy.  They  are  the  primary  standard 
candles  for  measuring  extragalactic  distances,  and  they  provide 
critical  tests  of  stellar-evolution  theory.  Surprisingly,  however, 
until  now  there  has  not  been  a  single  Cepheid  with  a  purely 
dynamical  measurement  of  its  mass. 

Polaris  ( a  UMi)  is  the  nearest  and,  at  second  magnitude, 
the  brightest  classical  Cepheid,  albeit  one  with  a  small  light 
amplitude  in  its  3.97-day  pulsation  period  (Turner  et  al.  2005  and 
references  therein).  The  amplitude,  which  had  been  decreasing 
for  several  decades,  now  appears  to  have  stabilized  and  may  be 
increasing  again  (Bruntt  et  al.  2008  and  references  therein).  The 
Hipparcos  parallax  of  Polaris  indicates  a  luminosity  consistent 
with  pulsation  in  the  first  overtone  (FO;  Feast  &  Catchpole 
1997). 

Polaris  is  the  brightest  member  of  a  triple  system  (see 
Kamper  1996  and  references  therein).  The  well-known  visual 
companion,  Polaris  B,  is  an  eighth  magnitude  F3  V  star  at  a 
separation  of  18".  The  Cepheid  itself  is  a  member  of  a  single- 
lined  spectroscopic  binary  with  a  period  of  30  yr.  In  this  paper, 
we  report  the  first  direct  detection  of  the  close  companion,  from 
which  we  derive  the  first  entirely  dynamical  mass  measurement 
for  a  Cepheid. 

Cepheid  masses  are  a  key  parameter  for  testing  stellar  evo¬ 
lutionary  calculations.  Beginning  in  the  1960s,  discrepancies 


*  Based  on  observations  with  the  NASA/ESA  Hubble  Space  Telescope 
obtained  at  the  Space  Telescope  Science  Institute,  which  is  operated  by  the 
Association  of  Universities  for  Research  in  Astronomy,  Inc.,  under  NASA 
contract  NAS5-26555. 


were  found  in  the  sense  that  Cepheid  masses  derived  from  pulsa¬ 
tion  modeling  were  lower  than  those  derived  from  evolutionary 
tracks.  A  revision  in  envelope  opacities  brought  pulsation  and 
evolutionary  masses  closer  together,  partially  alleviating  this 
“Cepheid  mass  problem.”  However,  recent  evolutionary  and 
pulsation  constraints  for  Galactic  (Bono  et  al.  2001b;  Caputo 
et  al.  2005;  Keller  2008;  Natale  et  al.  2008)  and  Magellanic 
(Bono  et  al.  2002;  Keller  &  Wood  2006)  Cepheids  still  imply 
a  discrepancy  in  masses  at  the  ~  15-20%  level.  The  luminosi¬ 
ties  of  Cepheids  depend  on  the  mass  of  the  helium-burning 
core;  physical  mechanisms  affecting  the  helium  core  mass  in¬ 
clude  mixing  due  to  convective  core  overshoot  during  the  main- 
sequence  phase,  mass  loss,  stellar  rotation,  and  radiative  opacity. 
A  directly  measured  mass  for  Polaris  would  thus  provide  an  im¬ 
portant  constraint  on  this  theoretical  framework. 

2.  OBSERVATIONS  AND  DATA  REDUCTION 

With  the  intention  of  a  direct  detection  of  the  close  com¬ 
panion,  we  imaged  Polaris  with  the  Hubble  Space  Telescope 
( HST)  and  the  High  Resolution  Channel  (HRC;  plate  scale 
0"026  pixel-1)  of  the  Advanced  Camera  for  Surveys  (ACS). 
We  chose  the  ultraviolet  (UV)  F220W  filter  (effective  wave¬ 
length  ~  2255  A)  in  order  to  minimize  the  contrast  between 
Polaris  and  the  close  companion,  which  we  anticipated  to  be  a 
main- sequence  star  slightly  hotter  than  the  Cepheid,  and  also  to 
minimize  the  size  of  the  point- spread  function  (PSF). 

Observations  were  obtained  on  2005  August  2-3  and  again 
on  2006  August  13.  At  the  first  epoch,  we  obtained  a  series  of 
0. 1-0.3  s  exposures  dithered  across  200  pixels  on  the  chip  over 
the  course  of  one  HST  orbit,  with  several  exposures  taken  at  each 
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Figure  1.  Co-added  ACS  HRC  images  of  Polaris  Aa,  Ab  taken  with  the  F220W  filter  on  2005  August  2  (left)  and  2006  August  13  (middle).  The  close  companion 
Ab  is  detected  at  the  lower  left  of  the  primary  (at  about  a  “7  o’clock”  position).  The  images  are  0785  x  0785  and  the  directions  of  north  and  east  are  indicated.  The 
right-hand  panel  shows  a  co-added  image  of  Polaris  B  from  longer  exposures  taken  during  the  2006  observations,  and  scaled  to  the  flux  level  of  the  Polaris  Aa,  Ab 
images.  There  is  no  artifact  in  the  Polaris  B  PSF  at  the  location  of  Ab. 
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Figure  2.  Contour  maps  of  the  co-added  images  shown  in  Figure  1.  The  outermost  contour  intervals  in  each  panel  range  from  0.001  to  0.010  of  the  peak  flux  in  steps 
of  0.001,  and  thereafter  are  at  levels  of  0.02,  0.04,  0.08,  0.16,  0.32,  and  0.64  of  the  peak  flux.  The  contours  again  demonstrate  the  absence  of  any  artifact  at  the  location 
of  the  Ab  companion. 


Figure  3.  Archival  images  of  flux-calibration  standard  stars  observed  with  the  ACS/HRC  in  the  F220W  filter.  The  star  names  and  dates  of  observation  are  listed  in 
each  panel.  There  is  no  artifact  at  the  location  of  the  Polaris  Ab  companion. 


dither  position.  At  the  second  epoch,  we  used  the  same  dither 
pattern,  but  divided  the  spacecraft  orbit  between  a  series  of  0.3  s 
exposures  on  Polaris  A  and  20  s  exposures  on  Polaris  B.  For  the 
longer  exposures,  we  placed  Polaris  B  at  the  same  chip  location 
as  Polaris  A  in  the  short  exposures,  so  as  to  provide  an  accurate 
PSF  for  a  single  star  at  the  same  place  in  the  field  of  view. 

Figure  1  shows  the  co-added  images  of  Polaris  A  from  2005 
(left  panel)  and  2006  (middle  panel).  The  close  companion 
(which  we  designate  Polaris  Ab)  is  detected  at  the  lower  left 
of  the  primary  (at  about  a  “7  o’clock”  position).  Because  of  the 
asymmetric  PSF  shape,  we  performed  several  checks  to  confirm 
that  the  apparent  companion  is  not  an  artifact.  The  right-hand 
frame  in  Figure  1  shows  Polaris  B  in  the  2006  image,  with  the 
star  shifted  to  the  same  field  position  as  Polaris  A,  and  with  its 


image  scaled  to  the  same  flux  as  Polaris  A.  This  PSF  shows  no 
artifact  at  the  location  of  the  companion  seen  in  the  images  of 
Polaris  A. 

Figure  2  shows  contour  maps  of  the  same  three  images.  Again 
the  faint  companion  is  seen  in  the  2005  and  2006  images  of 
Polaris  A,  and  there  is  no  PSF  artifact  at  this  location  in  the 
image  of  Polaris  B.  We  also  retrieved  images  from  the  HST 
archive  of  several  standard  stars  observed  with  ACS/HRC  in 
the  F220W  filter  over  an  interval  of  4  yr.  Examples  of  these 
observations  are  shown  in  Figure  3.  Although  the  PSF  structure 
does  vary  somewhat  with  time  due  to  small  changes  in  telescope 
focus  and  other  instrumental  phenomena  none  of  these  images 
shows  any  artifact  at  the  location  of  the  Polaris  Ab  companion 
seen  in  Figures  1  and  2. 
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Table  1 

Position  and  Magnitude  of  Polaris  Ab  Relative  to  Polaris  Aa 


Besselian  Date 

UT  Date  &  Time 

P(") 

PA.  (J2000)  (°)a 

Am(F220W) 

2005.5880 

2006.6172 

2005  Aug  2  23:45 

2006  Aug  16  22:01 

0.172  ±0.002 

0.170  ±0.003 

231.4  ±0.7 

226.4  ±  1.0 

5.38  ±  0.09 

5.40  ±  0.09 

Note.  a  RA.s  for  equinox  of  date  are  232?8  and  228°  1. 

Table  2 

Position  and  Magnitude  of  Polaris  B  Relative  to  Polaris  Aa 

Besselian  Date 

UT  Date  &  Time 

P  (") 

P.A.  (J2000)  (°)a 

Am(F220W) 

2005.5880 

2006.6172 

2005  Aug  2  23:45 

2006  Aug  16  22:01 

18.217  ±0.003 

18.214  ±0.003 

230.540  ±  0.009 

230.520  ±  0.009 

4.53  ±  0.04 

4.45  ±  0.02 

Note.  a  PA.s  for  equinox  of  date  are  231°980  and  232^216. 


Table  3 

Previously  Determined  Orbital  Parameters 


Kamper  (1996) 

Wielen  et  al.  (2000) 

Wielen  et  al.  (2000) 

Prograde 

Retrograde 

P  (yr) 

29.59  ±  0.02 

50.1  ±4.8 

in 

130.2  ±  4.9 

T 

e 

con 

KAa  (km  s-1) 

1987.66  ±  0.13 

0.608  ±  0.005 

303.01  ±  0.75 

3.72  ±  0.03 

Qa(°) 

276.2  ±  9.5 

Qa(°) 

167.1  ±9.4 

Note.  a  Values  of  Q  quoted  in  Wielen  et  al.  (2000)  correspond  to  the  astrometric  orbit  of  Polaris  Aa 
relative  to  the  center  of  mass. 


To  measure  the  separation  and  position  angle  between 
Polaris  Aa  and  the  newly  revealed  close  companion,  we  used 
the  calibrated  flat-fielded  exposures  provided  by  the  HST  re¬ 
duction  pipeline.  At  each  dither  location,  we  median- filtered 
the  repeated  observations  to  remove  cosmic  rays.  We  then  ex¬ 
tracted  subarrays  from  the  images  centered  on  Polaris  Aa  with 
a  size  of  0"85  x  0"85.  We  used  the  observations  of  Polaris  B 
from  2006  as  a  reference  PSF  to  construct  models  of  the  close 
pair  (Aa,  Ab)  by  searching  through  a  grid  of  separations  and 
flux  ratios.  The  IDL  interpolate  procedure  was  used  to  shift 
the  PSF  by  subpixel  intervals,  using  cubic  convolution  inter¬ 
polation.  The  background  was  least- squares  fitted  with  a  tilted 
plane. 

The  separations,  position  angles  (PA.s),  and  magnitude  dif¬ 
ferences,  determined  through  x2  minimization  between  the 
models  and  the  observations,  are  given  in  Table  1 .  The  uncertain¬ 
ties  were  determined  by  analyzing  multiple  images  individually 
and  computing  the  standard  deviation.  We  applied  the  filter- 
dependent  geometric  distortion  correction  of  Anderson  &  King 
(2004)  to  convert  the  pixel  values  to  a  separation  in  arcseconds. 
To  define  the  orientation  of  the  detector  y-axis  on  the  sky,  and 
thus  determine  the  PA.  of  the  binary  relative  to  the  pole  for 
the  equinox  J2000.0,  we  used  the  HST  image-header  keyword 
PA_APER. 

We  also  measured  the  separation  and  PA.  of  the  wide  com¬ 
panion,  Polaris  B,  relative  to  Aa,  and  the  results  are  presented 
in  Table  2.  The  good  agreement  with  the  historical  measure¬ 
ments  of  the  PA.  of  Polaris  B  relative  to  A  (see  Section  4) 
indicates  that  we  are  properly  defining  the  direction  of  north  in 
spite  of  the  extreme  northerly  declination.  (Since  the  historical 
double-star  convention  is  to  give  the  PA.  for  the  equinox  of  the 
date  of  observation,  we  computed  the  precession  corrections 
and  gave  the  adjusted  PA.s  in  the  footnotes  to  Tables  1  and  2, 
for  the  convenience  of  archivists.) 


3.  ORBITAL  SOLUTIONS 

We  stress  that  the  orbital  analyses  discussed  below  are  only 
preliminary  fits  to  data  with  a  very  limited  sample  (only  two 
points)  of  separations  and  PA.s.  We  followed  three  different 
approaches  to  determining  the  orbital  parameters,  in  order  to 
illustrate  the  scope  of  the  available  data. 

Kamper  (1996)  rederived  the  single-lined  spectroscopic  orbit 
of  Polaris  Aa  using  improved  radial-velocity  data,  and  a  careful 
removal  of  the  velocity  signal  due  to  the  Cepheid  pulsation.  His 
solution  provides  the  period,  time  of  periastron  passage,  eccen¬ 
tricity,  angle  between  the  node  and  periastron,  and  the  radial- 
velocity  semiamplitude  of  the  primary  (denoted  P,  T ,  e ,  co, 
and  ^Aa,  respectively).  By  comparing  the  Hipparcos  proper 
motion  of  Polaris  Aa  (which,  over  the  duration  of  the  Hipparcos 
mission,  is  nearly  instantaneous  in  the  context  of  the  30-year 
orbit)  with  the  ground-based  long-term  average  proper  motion 
from  the  FK5  (which  is  essentially  the  motion  of  the  center  of 
mass),  Wielen  et  al.  (2000)  determined  the  inclination  and  the 
PA.  of  the  line  of  nodes  (i  and  Q).  Their  analysis,  however, 
allows  for  retrograde  and  prograde  orbital  solutions  (the  two  or¬ 
bits  being  tangential  at  the  Hipparcos  epoch),  with  two  different 
values  of  i  and  £2.  The  orbital  parameters  based  on  the  Kamper 
(1996)  and  the  two  Wielen  et  al.  (2000)  solutions  are  presented 
in  Table  3. 

The  HST  detection  of  the  close  companion  Ab  and  its  orbital 
motion  at  two  epochs  establishes  a  retrograde  sense  for  the  orbit 
(thus  confirming  the  strong  preference  stated  by  Wielen  et  al.  for 
their  retrograde  solution).  Additionally,  it  provides  constraints 
on  the  remaining  unknown  parameter  of  the  orbit,  the  semimajor 
axis  a.  A  combination  of  the  spectroscopic  mass  function, 

f(M)  =  (MAb  sin/)3 /(Mm  +  MAb)2 

=  3.784  x  l(T5/s±,P(l  -  e2)3/2  , 
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Table  4 

Orbital  Parameters*1  and  Massesb 


Parameter 

Wielen  et  al.  (2000) 

Tit  HST 

Joint  Fit  to  HST 

Retrograde0 

Only 

and  Proper  Motion 

in 

130.2  (F) 

155^46 

128+^\ 

Qd  (°) 

347.1  (F) 

9-OtS 

19+1/ 

o 

0.131  ±  0.04 

0-1 16-oX 

0.133!™*5! 

Mtot(M0) 

5.6  ±  1.0 

O  Q+1.0 

MaAMq) 

4.3  ±  1.3 

2  1+L4 

Z.l— 2  i 

4  5+2-2 

MaUMq) 

1.26  ±  0.80 

1  8+1-8 
L8-0.5 

1  26+014 
i'Z0— 0.07 

Notes. 

a  Spectroscopic  orbital  parameters  (P,  T,  e,oo )  held  fixed  at  Kamper  (1996) 
values. 

b  Assumes  Hipparcos  parallax  of  7.72  ±  0.12  mas. 

c  Orbital  parameters  marked  (F)  are  fixed  at  the  given  values  when  computing 
the  best  fit  solution. 

d  Q  has  been  rotated  by  180°  from  the  values  quoted  in  Wielen  et  al.  (2000)  to 
correspond  to  the  orbit  of  Polaris  Ab  relative  to  Polaris  Aa. 


with  the  total  mass  from  Kepler’s  Third  Law, 

cp* 

Mot  =  M Aa  +  M Ab  =  3  2  , 

7 riPz 

then  yields  the  masses  of  the  binary  components, 


MAh  =  0.03357 


KAdia Vl  -  e2 

jz2P  sin  i 


(1) 


M\a  =  Mot  —  Mvb, 

where  a  and  the  parallax  n  are  in  arcseconds,  P  is  in  years,  KAsl 
is  in  km  s-1,  and  the  masses  are  in  M©. 

In  Sections  3. 1-3.3  we  describe  the  orbital  fits  that  we  com¬ 
puted  based  on  a  synthesis  of  the  spectroscopic,  astrometric,  and 
HST  data.  In  these  sections,  we  present  three  successively  more 
comprehensive  orbital  solutions.  Table  4  summarizes  the  orbital 
parameters  determined  from  each  of  these  fits;  the  individual 
columns  are  described  in  more  detail  in  Sections  3. 1-3.3. 

3.1.  Semimajor  Axis 

As  a  first  approximation  to  an  orbital  solution,  we  fixed 
the  spectroscopic  and  astrometric  parameters  ( P,T,e,co,i,Cl ) 
to  be  those  determined  by  Kamper  (1996)  and  Wielen  et  al. 
(2000),  and  solved  only  for  the  semimajor  axis  a  based  on 
the  two  HST  separation  measurements.  The  orbital  parameters 
for  this  solution  are  listed  in  the  second  column  of  Table  4. 
Figure  4  compares  this  retrograde  orbit  fit  with  the  HST  mea¬ 
surements,  and  shows  extremely  poor  agreement  for  the  P.A.s. 
However,  the  relatively  large  uncertainties  in  i  and  Q  deter¬ 
mined  by  Wielen  et  al.  (2000)  provide  considerable  flexibility 
for  adjusting  the  orbital  parameters  in  order  to  improve  the  fit 
quality. 


3.2.  Best  Fit  to  the  HST  Measurements 

To  get  a  better  fit  to  the  HST  data,  we  solved  for  i,  Q,  and 
a  based  on  the  two  separation  and  P.  A.  measurements  from  the 
HST  observations,  while  holding  the  relatively  well-determined 
spectroscopic  parameters  ( P,T,e ,  and  co)  fixed.  We  computed 
the  orbit  fit  through  a  standard  Newton-Raphson  method  in 
X2  space,  and  present  the  results  in  the  third  column  of  Table  4. 


0.1  0.0  -0.1  -0.2 
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Figure  4.  Comparison  of  the  retrograde  orbital  solution  of  Wielen  et  al.  (2000; 
solid  ellipse)  with  our  HST  measurements  of  Polaris  Ab  (filled  circles).  The 
arrow  indicates  the  direction  of  motion.  The  orbit  was  calculated  by  fixing 
the  spectroscopic  and  astrometric  orbital  parameters  (see  Table  3)  and  solving 
only  for  the  semimajor  axis  based  on  the  HST  separation  measurements.  We 
found  a  =  07'  131.  The  dotted  lines  connect  the  observed  positions  to  x  symbols 
marking  the  predicted  positions,  and  show  the  inadequacy  of  this  simple  solution 
at  predicting  the  P.A.s.  The  size  of  the  HST  error  bar  is  indicated  by  the  cross  at 
the  top  of  the  plot. 


To  explore  the  range  of  orbital  parameters  that  fit  the  HST 
data,  we  performed  a  Monte  Carlo  search  by  selecting  values 
of  i,  Cl,  and  a  at  random.  We  searched  for  10,000  solutions 
within  the  3a  confidence  interval,  corresponding  to  a  difference 
of  A/2  =  9  from  the  minimum  x2  value.  Figure  5  shows  cross¬ 
cuts  through  the  x2  surfaces  for  the  three  derived  parameters. 
Using  the  recently  revised  Hipparcos  parallax  of7.72±0.12 
mas  (van  Leeuwen  et  al.  2007),  we  computed  the  total  mass  of 
the  binary  through  Kepler’s  Third  Law  for  all  of  the  solutions 
found  in  the  Monte  Carlo  search.  In  the  last  panel  of  Figure  5,  we 
show  a  plot  of  the  total  mass  versus  inclination.  Because  a  visual 
orbit  is  insensitive  to  the  individual  masses  of  the  components, 
when  combining  the  total  mass  with  the  spectroscopic  mass 
function,  there  exist  values  of  the  inclination  that  produce 
negative  masses  for  MA3l  or  MA b.  In  the  remaining  analysis, 
we  removed  these  negative-mass  solutions  from  our  sample 
of  possible  orbits.  Essentially,  this  rejects  all  orbital  solutions 
with  i  >  168°.  The  la  uncertainties  listed  in  the  third  column 
of  Table  4  are  determined  from  the  Ax2  =  1  confidence 
interval  of  the  modified  distribution.  The  values  of  i  and  Q 
agree  with  the  retrograde  parameters  computed  by  Wielen 
et  al.  (2000)  at  the  1.5a  and  2.2a  levels  respectively.  Figure  6 
shows  three  examples  illustrating  how  the  orbit  fit  varies  within  a 
la  confidence  interval. 

We  note  that  the  HST- only  solution  yields  a  secondary 
mass  of  MAb  =  1.8  AT©  (Table  4,  column  3),  corresponding 
approximately  to  an  A5  V  star.  The  discussion  in  Section  5.1,  as 
well  as  the  lack  of  a  detection  of  the  companion  in  UV  spectra 
obtained  by  Evans  (1988)  with  the  International  Ultraviolet 
Explorer  ( IUE ),  makes  it  highly  improbable  that  Ab  could  be 
this  hot.  With  only  two  HST  measurements  sampling  the  orbit 
thus  far,  we  do  not  yet  have  a  good  constraint  on  the  curvature, 
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Figure  5.  Crosscuts  through  the  x2  surfaces  derived  from  the  fit  to  the  HST  measurements  alone  while  holding  the  spectroscopic  parameters  fixed.  These  results  were 
obtained  by  performing  a  Monte  Carlo  search  for  orbital  solutions  within  the  3 o  confidence  interval  (Ax2  =  9).  The  total  mass  was  derived  for  each  of  the  10,000 
orbits  found  within  this  interval  using  a  parallax  of  7.72  mas.  The  dotted  lines  indicate  the  critical  value  of  the  inclination  (i  =  168°),  above  which  negative  values  of 
a  component  mass  are  produced.  The  color  codes  in  the  electronic  version  of  this  figure  correspond  to  the  1  o  (red),  2cr  (blue),  and  3cr  (green)  confidence  intervals. 


0.1  0.0  -0.1  -0.2 
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Figure  6.  Three  examples  of  orbits  resulting  from  variations  of  the  orbital 
parameters  within  the  lcr  confidence  intervals  based  on  the  fit  to  the  HST 
measurements  alone  while  holding  the  spectroscopic  parameters  fixed.  The 
arrows  indicate  the  direction  of  motion.  All  three  orbits  fit  the  HST  measurements 
(filled  circles)  within  the  lcr  error  bar  shown  at  the  top  of  the  figure,  but  they 
imply  total  system  masses  ranging  from  3.6  to  4.9  Af0. 

and  hence  acceleration,  of  the  orbit.  In  turn,  this  limits  how  well 
we  can  determine  the  total  system  mass,  and  it  contributes  to 
the  large  errors  quoted  in  Table  4. 


3.3.  Joint  Fit  to  HST  and  Proper-Motion  Measurements 

Incorporation  of  the  Hipparcos  proper-motion  measurements 
into  the  orbital  fit  extends  the  time  coverage  of  the  measurements 
to  ~15  yr.  This  represents  a  significant  fraction  of  the  orbital 
period  and  is  therefore  likely  to  improve  the  reliability  of  the 
results.  Following  the  technique  described  in  Wielen  et  al. 
(2000),  we  performed  a  simultaneous  fit  to  the  proper-motion 
data  and  the  HST  measurements. 

As  Wielen  et  al.  point  out,  the  FK5  proper  motion  is  aver¬ 
aged  over  several  cycles  of  the  orbital  period  and  therefore  re¬ 
flects  the  center-of-mass  motion  of  Polaris  Aa,  Ab.  Because  of 
the  shorter  time  span  of  the  Hipparcos  mission,  the  Hipparcos 
proper  motion  more  nearly  represents  an  instantaneous  mea¬ 
surement  of  the  combined  proper  motion  of  the  center  of  mass 
of  the  Aa,  Ab  pair  and  the  orbital  motion  of  the  photocenter  about 
the  center  of  mass  at  the  epoch  of  the  observations  (~  199 1.3). 
The  difference  between  the  FK5  and  Hipparcos  proper  motions 
thus  gives  the  offset  caused  by  the  orbital  motion. 

In  computing  the  joint  fit,  we  held  the  spectroscopic  param¬ 
eters  (P,  T,  e,  co,  and  Ka&)  fixed  at  the  Kamper  (1996)  values 
and  solved  for  i,  Q,  and  a,  again  using  a  Newton-Raphson 
technique  in  x2  space.  The  input  data  were  the  relative  posi¬ 
tions  of  Aa  and  Ab  at  the  two  HST  epochs,  and  the  difference 
between  the  Hipparcos  and  FK5  proper  motions.  To  incorpo¬ 
rate  the  proper-motion  data  into  the  orbit  fit,  we  had  to  com¬ 
pute  the  time-dependent  offset  of  the  photocenter  relative  to  the 
center  of  mass  predicted  by  the  orbital  parameters  during  the 
time  of  the  Hipparcos  observations.  To  compute  these  offsets, 
we  converted  the  semimajor  axis  of  Polaris  Aa  determined  from 
the  single-lined  spectroscopic  orbit  to  the  semimajor  axis  of  the 
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Figure  7.  Schematic  comparison  between  the  instantaneous  and  mean  proper 
motions  of  Polaris  A.  The  vectors  indicate  the  magnitude  and  direction  of 
the  annual  proper-motion  measurements  from  Hipparcos  (/xhip;  red),  the 
long-term  average  FK5  ground-based  motion  (/xfks;  blue),  and  the  vector 
difference  between  the  two  (A /x;  black).  Dotted  red  and  blue  lines  indicate 
the  lcr  uncertainties  in  the  Hipparcos  and  FK5  measurements.  The  green  vector 
represents  the  best  fit  difference  between  the  proper  motions  (A/Xfit)  computed 
from  our  simultaneous  orbit  fit  to  the  HST  measurements  and  the  proper-motion 
data  while  holding  the  spectroscopic  parameters  fixed  (see  Section  3.3). 


Table  5 

Proper-Motion  Data  for  Polaris  A 


Quantity 

System 

/xa  cos  S 
(mas  yr-1) 

/x<$ 

(mas  yr-1) 

T  FK5 

FK5 

+38.30  ±  0.23 

-15.20  ±  0.35 

TFK5 

HIP 

+41.50  ±  0.97 

-16.73  ±  0.75 

/XHIP 

HIP 

+44.22  ±  0.47 

-11.74  ±  0.55 

A/x  =  /X HIP  -  /x FK5 

HIP 

+2.72  ±  1.08 

+4.99  ±  0.93 

A/Xfit 

HIP 

+4.59  ±  2.52 

+1.21  db  0.74 

photocenter  by  using  the  mass  ratio  computed  from  the  full  set  of 
orbital  parameters  and  a  magnitude  difference  between  Aa  and 
Ab  of  AV  =  7.2  (see  Section  5.1).  This  conversion  is  specified 
by  Equations  (8)  and  (9)  of  Wielen  et  al.  (2000).  The  expected 
difference  between  the  instantaneous  and  the  mean  proper  mo¬ 
tion  (A/Xfit)  at  the  central  epoch  of  the  Hipparcos  observations 
is  then  computed  from  Equation  (18)  of  Wielen  et  al. 

Table  5  shows  the  values  of  the  proper  motions  used  by 
Wielen  et  al.  (2000).  The  first  line  shows  the  proper  motion 
of  Polaris  given  in  the  FK5  catalog  (/xfks)-  In  the  second  line, 
a  systematic  correction  is  applied  to  convert  the  proper  motion 
from  the  FK5  reference  system  to  the  Hipparcos /International 
Celestial  Reference  System  (ICRS)  system  (Wielen  et  al.  2000). 
The  proper  motion  measured  by  Hipparcos  (/xhip)  is  given  in 
the  third  line.  The  difference  in  the  proper  motions  measured 
A/x  =  /x hip  —  /xfk5  is  given  in  the  fourth  line. 

The  last  line  in  Table  5  shows  the  best-fit  difference  between 
the  instantaneous  and  mean  proper  motions  (A/Xfit)  calculated 
from  our  simultaneous  fit  to  the  relative  separation  and  P.A. 
measurements  of  Polaris  Aa,  Ab  and  the  proper-motion  data. 
Figure  7  shows  a  graphical  representation  of  the  Hipparcos  and 
FK5  proper  motions,  and  compares  our  best-fit  value  with  the 
measured  value  for  the  difference  between  the  instantaneous 
and  mean  proper  motions. 

As  Figure  7  illustrates,  the  best-fit  value  of  A/x  in  right  as¬ 
cension  agrees  well  (0.7 a)  with  the  measured  value,  but  the 
agreement  is  poorer  (3.2a)  in  declination.  This  discrepancy 
probably  arises  from  our  constraining  the  spectroscopic  param¬ 
eters  to  be  exactly  those  derived  by  Kamper  (1996),  thus  forcing 
the  proper  motions  and  HST  measurements  to  absorb  the  errors. 
We  found  that  by  allowing  some  variation  in  the  spectroscopic 
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Figure  8.  Final  best  fit  orbit  of  Polaris  Ab  relative  to  Aa  based  on  a  simultaneous 
fit  to  the  HST  measurements  (filled  circles)  and  the  proper-motion  data  while 
holding  the  spectroscopic  parameters  fixed.  The  predicted  orbital  positions  at 
the  epochs  of  the  HST  measurements  are  marked  by  the  x  symbols  (partially 
hidden  by  the  observed  points).  The  position  of  the  companion  during  the  time 
frame  of  the  Hipparcos  mission  is  highlighted  by  the  shaded  gray  line  segment; 
its  direction  of  motion  is  of  course  180°  different  from  the  direction  of  A/Xfit, 
shown  as  a  green  arrow  in  Figure  7. 

parameters,  we  could  substantially  improve  the  fit  to  the  proper- 
motion  measurements.  Once  we  have  sampled  enough  of  the 
visual  orbit  to  better  constrain  /,  Q,  and  a,  the  optimal  orbital 
solution  should  be  found  by  doing  a  simultaneous  fit  to  the  radial 
velocity,  HST ,  and  proper-motion  data.  Unfortunately,  however, 
Kamper  only  tabulated  the  radial  velocities  before  removal  of 
the  pulsational  variation,  so  a  re-computation  of  the  pulsation 
corrections  would  have  to  be  carried  out — a  task  well  beyond 
the  scope  of  this  paper  and  one  that  should  await  the  availability 
of  more  HST  observations. 

The  orbital  parameters  and  derived  masses  from  our  final 
combined  fit  are  given  in  the  last  column  of  Table  4,  which 
also  contains  our  final  best  estimates  for  the  dynamical  masses 
of  both  stars.  The  best-fit  orbit  of  Polaris  Ab  relative  to 
Aa  is  plotted  in  Figure  8.  The  gray-shaded  portion  of  the  orbit 
marks  the  location  of  the  companion  during  the  interval  of  the 
Hipparcos  observations  in  the  early  1990s,  and  it  should  be 
noted  that  the  direction  of  motion  at  that  time  was,  of  course, 
180°  different  from  the  direction  of  the  differential  motion  of 
Polaris  Aa  (A/Xfit)  shown  as  a  green  arrow  in  Figure  7. 

4.  ASTROMETRY  OF  POFARIS  B 

Visual  measurements  of  the  P.A.  and  separation  of  the  wide 
companion  Polaris  B  relative  to  Polaris  A  extend  back  to  the 
early  nineteenth  century,  with  a  few  photographic  observations 
that  are  available  from  the  twentieth  century.  We  have  compiled 
these  measurements  from  Kamper  (1996)  and  the  Washington 
Double  Star  Catalog  (Mason  et  al.  2001).7 

In  Figure  9  we  plot  the  P.A.  measurements  of  Polaris  B,  which 
have  been  precessed  to  the  J2000.0  equinox.  Due  to  its  slow 
relative  motion  and  large  magnitude  difference,  Polaris  A,  B 


7  http://ad.usno.navy.mil/wds 
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Figure  9.  Historical  measurements  of  the  position  angle  of  Polaris  B  relative 
to  Polaris  A,  precessed  to  the  equinox  of  J2000.0.  Our  HST  measurements  are 
marked  by  the  two  filled  circles.  The  solid  line  shows  a  least-squares  fit  yielding 
a  rate  of  change  in  the  position  angle  of  —  0°00035  ±  0?00094  yr-1 . 


was  generally  not  included  in  the  major  double-star  observing 
programs,  especially  in  the  latter  half  of  the  twentieth  century. 
The  result  of  this  was  that  most  of  the  measures  that  were 
obtained  tended  to  be  made  by  less-experienced  observers, 
often  using  older  techniques  and/or  smaller  telescopes;  this 
may  explain  the  surprisingly  large  scatter  in  the  late  twentieth 
century.  A  linear  least- squares  fit  to  the  data  yields  a  rate  of 
change  in  the  PA.  of  — (F00035  =b  (T00094  yr-1,  consistent 
with  no  detectable  change  in  PA.  for  the  past  two  centuries. 
(The  earlier  work  of  Kamper  had  given  a  marginal  detection  of 
+0?0086  zb  0?0076  yr-1.) 

The  separation  measurements  for  Polaris  B  are  plotted  in 
Figure  10.  There  has  been  a  slow  downward  trend  in  the 
separation,  with  a  least-squares  fit  giving  a  rate  of  —  1 .67  ±  0. 19 
mas  yr-1.  From  ground-based  measures  only,  Kamper  (1996) 
had  found  —1.7  zb  0.6  mas  yr-1.  Since  the  absolute  proper 
motion  of  Polaris  A  is  ~46  mas  yr-1  (see  Table  5),  the  absolute 
motions  of  A  and  B  agree  to  within  about  4%.  At  the  distance 
of  Polaris  given  by  the  Hipparcos  parallax,  the  difference  in 
tangential  velocities  between  A  and  B  is  1.03  =b  0.12  km  s-1. 

In  computing  these  least-squares  fits,  we  weighted  the  ob¬ 
servations  by  estimates  of  their  measurement  errors.  For  the 
HST  observations,  we  applied  the  uncertainties  quoted  in 
Table  2.  For  the  historical  measurements,  we  divided  the  data 
into  four  groups  spanning  approximately  50  years  each.  We 
assumed  measurement  uncertainties  equal  to  the  standard  devi¬ 
ation  of  the  values  measured  in  each  of  these  four  groups. 

The  slowly  diminishing  separation  of  A  and  B  (at  constant 
PA.)  is  not  inconsistent  with  orbital  motion  in  a  physically 
bound  pair — which  is  also  supported  by  the  close  similarities 
of  the  radial  velocities  of  A  and  B  (Kamper  1996;  Usenko 
&  Klochkova  2008).  To  predict  an  order-of-magnitude  rate  of 
change  in  the  separation,  we  assumed  a  circular  orbit  with  a 
total  system  mass  of  MAa  +  MAb  +  MB  =7.15  M©  (based  on 
the  Aa+Ab  mass  of  5.8  M©  in  the  last  column  of  Table  4  and 
a  mass  for  Polaris  B  of  1.35  M© — see  below).  Adopting  the 
revised  Hipparcos  parallax  oil .12  mas,  and  assuming  an  edge- 
on  orbit  with  a  period  of  ~  100,000  yr,  we  find  a  semimajor  axis 
of  a  ~  32"  (or  0.02  pc).  At  the  orbital  phase  implied  by  the 
observed  separation  of  1 8/2,  the  relative  motion  would  then  be 
—  1.65  mas  yr-1,  close  to  the  observed  value. 


Figure  10.  Separation  measurements  of  Polaris  B  relative  to  Polaris  A. 
Our  HST  measurements  are  marked  by  the  two  filled  circles.  The  solid  line 
shows  a  linear  least-squares  fit  yielding  a  rate  of  change  in  the  separation  of 
— 1.67  ±  0.19  mas  yr-1. 


Table  6 

Spectral-Type  Comparisons 


Star 

Spectral 

type 

Mass 

(Mo) 

Parallax 

(mas) 

Mv 

F220W  flux 

relative  to  Polarisa 

78  UMa 

F2  V 

1.41 

40.06 

+2.9 

0.0262 

Polaris  B 

0.0178 

HD  27524 

F5  V 

1.33 

19.55 

+3.2 

0.0114 

Polaris  Ab 

0.0074 

HD  27808 

F8  V 

1.22 

24.47 

+4.1 

0.0029 

Note.  aFlux  ratios  for  78  UMa,  HD  27524,  and  HD  27808  are  predicted 
from  their  IUE  spectra,  scaled  to  the  distance  of  Polaris;  ratios  for  Polaris  B 
and  Ab  are  those  observed  by  us,  corrected  for  red  leak  as  described  in  the 
text. 

5.  ASTROPHYSICAL  PROPERTIES  OF  THE 
COMPANIONS  OF  POLARIS 

5.7.  Polaris  Ab 

Our  observed  UV  magnitude  difference  between  Polaris  Aa 
and  Ab  may  be  used  to  infer  the  spectral  type,  and  hence  the 
mass,  of  the  newly  resolved  close  companion. 

We  downloaded  UV  spectra  from  the  IUE  data  archive8 
for  three  F-type  dwarfs  having  accurate  spectral  types  and 
parallaxes,  as  well  as  for  Polaris  itself.  The  F  stars  selected 
were  78  UMa  (HR  4931,  HD  1 13139;  F2  V),  HD  27524  (F5  V), 
and  HD  27808  (F8  V).  The  latter  two  stars  (Hyades  members), 
as  well  as  Polaris  itself,  were  taken  to  be  unreddened,  while  the 
spectrum  of  78  UMa  (a  member  of  the  Ursa  Major  group)  was 
dereddened  by  E(B  —  V)  =  0.01  mag.  We  then  scaled  the  flux 
distributions  for  the  three  stars  to  the  distance  of  Polaris,  using 
the  respective  Hipparcos  parallaxes. 

In  Table  6,  we  list  these  comparison  stars,  their  spectral 
types,  masses  implied  by  the  spectral  types,  their  Hipparcos 
parallaxes,  absolute  magnitudes  based  on  the  parallaxes,  and 
finally  the  predicted  flux  ratios  relative  to  Polaris  in  the 
ACS/HRC  F220W  band.  The  adopted  relationship  between 
spectral  types  and  masses  is  that  of  Harmanec  (1988).  The 


8  The  IUE  data  were  obtained  from  the  Multimission  Archive  at  the  Space 
Telescope  Science  Institute  (MAST).  Support  for  MAST  for  non -HST  data  is 
provided  by  the  NASA  Office  of  Space  Science  via  grant  NAG5-7584  and  by 
other  grants  and  contracts. 
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Figure  11.  IUE  UV  spectra  of  Polaris  and  of  three  main-sequence  F-type  stars  scaled  to  the  distance  of  Polaris  using  the  respective  Hipparcos  parallaxes.  The  two 
horizontal  lines  correspond  to  the  mean  flux  levels  in  the  F220W  filter  of  the  ACS/HRC  of  Polaris  B  (top)  and  Ab  (bottom).  The  lengths  of  the  horizontal  lines 
correspond  to  the  FWHM  of  the  filter.  On  the  basis  of  its  UV  flux  level,  Polaris  B  is  inferred  to  have  a  spectral  type  near  F3-F4  V,  in  good  agreement  with  its 
ground-based  classification  at  F3  V.  Polaris  Ab  is  inferred  from  its  F220W  flux  to  have  a  spectral  type  near  F6  V. 


flux  ratios  were  calculated  by  convolving  the  F220W  system- 
throughput  function  (see  Chiaberge  &  Sirianni  2007)  with  the 
scaled  IUE  spectra,  and  then  ratioing  with  respect  to  Polaris.  In 
Figure  1 1,  we  show  the  IUE  spectrum  of  Polaris  and  the  scaled 
spectra  of  the  three  F  dwarfs. 

As  listed  in  Table  1,  the  observed  magnitude  difference  in 
the  ACS/HRC  F220W  filter  between  Polaris  Aa  and  Ab  is 
5.39  =b  0.08  mag,  or  a  flux  ratio  of  0.0070  zb  0.0006.  A  small 
correction  to  this  ratio  is  needed  because  of  the  small  (~10- 
15%)  contribution  to  the  signal  from  the  red  leak  in  the  F220W 
filter,  Polaris  Aa  being  slightly  redder  than  the  companion. 
The  red-leak  contributions  have  been  tabulated  as  a  function 
of  spectral  type  by  Chiaberge  &  Sirianni  (2007),  leading  to 
a  corrected  in-band  flux  ratio  of  0.0074  zb  0.0006.  This  value 
is  entered  in  the  fourth  row  of  Table  6,  and  is  marked  with  a 
horizontal  line  in  Figure  1 1 .  Interpolation  in  the  last  column  of 
Table  6  then  leads  to  an  inferred  spectral  type  of  about  F6  V, 
an  absolute  visual  magnitude  of  Mv  —  +3.6,  and  an  expected 
mass  of  1.3  M0. 

The  apparent  V-band  magnitude  of  Polaris  Ab,  for  a  distance 
modulus  (m  —  M) o  =  5.56,  is  inferred  to  be  about  9.2,  or 
some  7.2  mag  fainter  at  V  than  Polaris  Aa.  This  illustrates  the 
advantage  of  observing  the  Polaris  system  in  the  UV,  which 
lessens  the  contrast  by  nearly  2  mag. 

5.2.  Polaris  B 

As  listed  in  Table  2,  we  also  measured  the  F220W  flux 
difference  between  Polaris  Aa  and  B  as  4.49  zb  0.04  mag,  or 
a  flux  ratio  of  0.0160  zb  0.0006.  Correction  for  red  leak,  as 
described  above,  changes  the  flux  ratio  to  0.0178  zb  0.0006, 
entered  in  the  second  row  of  Table  6,  and  also  marked  with  a 
horizontal  line  in  Figure  1 1 .  Interpolating  again  in  the  table,  we 
see  that  this  ratio  corresponds  to  a  star  intermediate  between 
types  F3  V  and  F4  V.  In  the  case  of  the  well-resolved  Polaris  B, 
the  optical  spectral  type  has  been  determined  from  the  ground. 
Our  result  is  in  gratifying  agreement  with  the  spectral  type  of 
F3  V  found  by  Turner  (1977)  and  Usenko  &  Klochkova  (2008), 
who  also  cite  earlier  spectral  classifications  of  Polaris  B  by 
experts  such  as  Bidelman.  This  finding  not  only  validates  our 


photometric  analysis  of  Polaris  Ab  above,  but  also  supports 
the  physical  association  of  Polaris  A  with  B.  Based  on  the 
relationship  between  the  spectral  type  and  mass  in  Table  6, 
we  infer  the  mass  of  Polaris  B  to  be  near  1.35  M©. 

Using  the  same  method  as  for  Polaris  Ab,  we  can  use  the 
UV  flux  ratio  to  infer  the  V  magnitude  of  Polaris  B  to  be  8.7. 
The  visual  magnitude  of  Polaris  B  can  be  measured  from  the 
ground,  but  is  made  difficult  by  scattered  light  from  Polaris  A. 
Kamper  (1996)  used  CCD  imaging  to  determine  a  magnitude 
difference  with  respect  to  A  of  A  V  =  6.61,  implying  V  =  8.59 
(in  good  agreement  with  earlier  photoelectric  measurements 
of  V  =  8.5  and  8.60  by  Fernie  1966  and  McNamara  1969, 
respectively;  Fernie  included  an  approximate  correction  for 
scattered  light  and  McNamara  states  that  he  observed  only  on 
excellent  nights).  In  more  recent  work,  to  be  reported  separately, 
we  have  been  carrying  out  astrometry  of  Polaris  B  with  the  Fine 
Guidance  Sensors  (FGS)  onboard  HST.  As  a  byproduct,  these 
observations  yield  an  accurate  V  magnitude  of  8.65  zb  0.02. 
Thus,  our  indirectly  inferred  V  magnitude  for  Polaris  B  of  8.7 
agrees  very  well  with  the  ground-  and  HST- based  observations. 

6.  DYNAMICAF  MASSES 
6.1.  Polaris  Ab 

The  final  column  in  Table  4  lists  the  dynamical  masses  of 
both  components  of  the  close  pair  Aa,  Ab  obtained  from  our 
final  orbital  solution,  as  described  in  Section  3.3.  For  Ab,  the 
dynamical  mass  is  1 .26^007  M©.  This  is  in  remarkably  good 
agreement  with  1.3  Mq  inferred  indirectly  from  the  UV  flux 
difference  (Section  5.1),  and  is  an  indicator  of  the  validity  of 
our  orbital  solution. 

6.2.  Theoretical  Implications  of  the  Cepheid’s  Dynamical  Mass 

The  dynamical  mass  of  the  Cepheid  Polaris  Aa  from  our 
final  orbital  solution,  as  listed  in  the  last  column  in  Table  4,  is 
4.5Z2124M0. 

We  compare  this  result  first  with  theoretical  “evolutionary” 
masses,  Me.  The  input  data  are  the  intensity- averaged  mean 
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Figure  12.  Evolutionary  tracks  (see  Pietrinferni  et  al.  2006)  in  the  My  vs. 
(B  —  V)o  plane,  with  the  location  of  Polaris  marked  (open  triangle  enclosing 
error  bar).  Top  panel:  canonical  evolutionary  models  neglecting  convective 
core  overshooting  during  H-burning  phases.  Bottom  panel:  noncanonical 
evolutionary  models  including  convective  core  overshooting.  The  assumed 
Polaris  distance  and  helium  (7)  and  metal  (Z)  abundances  (scaled  solar)  are 
indicated  in  the  figure,  and  the  legends  on  the  left  indicate  the  stellar  masses  for 
each  track.  These  models  include  mass  loss  with  a  Reimers  parameter  r]  =  0.4. 

apparent  magnitudes  {my  —  1.98,  mg  =  2.58,  from  Femie  et  al. 
1995  and  assumed  to  be  unreddened),  the  revised  Hipparcos 
distance  of  129.5  ±  2.0  pc  (van  Leeuwen  et  al.  2007),  and  a 
solar  metal  abundance  (Luck  &  Bond  1986;  Usenko  et  al.  2005). 
We  adopt  the  mass-period-luminosity  (MPL)  relation  for  He- 
burning  fundamental  pulsators  provided  by  Caputo  et  al.  (2005, 
their  Table  4).  Before  using  this  relation,  we  fundamentalized  the 
FO  pulsation  period  of  Polaris  with  the  relationship  log  Py  = 
logPpo  +  0.13.  By  assuming  Cepheid  luminosities  predicted 
by  “canonical”  evolutionary  models  that  neglect  convective- 
core  overshooting,  we  find  Me  =  6. 1  ±  0.4  M0.  However,  if  we 
assume  luminosities  predicted  by  “noncanonical”  evolutionary 
models  that  account  for  mild  convective-core  overshooting, 
given  by  L/Lcan  ~  1.3,  we  find  Me  =  5.6  ±  0.4  Af0.  Using 
the  mass-color-luminosity  (MCL)  relation  (Caputo  et  al.  2005, 
Table  5)  yields  very  similar  evolutionary  masses. 

The  H-R  diagrams  in  Figure  12  show  a  direct  comparison 
between  theoretical  evolutionary  tracks  and  observations  in  the 
My ,  ( B  —  V)o  plane.  In  both  panels  of  Figure  12,  we  plot 
the  location  of  Polaris  with  an  open  triangle  enclosing  a  small 
error  bar.  The  top  panel  shows  canonical  evolutionary  tracks  at 


solar  chemical  composition;  Me  ~  6  M0  provides  a  good  fit  to 
the  position  of  Polaris.  The  bottom  panel  shows  noncanonical 
tracks,  suggesting  Me  ~  5.5  M0,  except  that  the  tip  of  the 
blue  loop  is  not  quite  as  hot  as  Polaris.  However,  the  blueward 
extension  of  the  loops  is  affected  by  chemical  composition  and 
by  physical  and  numerical  assumptions  (Stothers  &  Chin  1991; 
Chiosi  et  al.  1992;  Bono  et  al.  2000;  Meynet  &  Maeder  2000; 
Xu  &  Li  2004). 

To  compare  our  result  with  “pulsation”  masses,  Mp,  we  used 
the  mass-dependent  period-luminosity-color  (PLC)  relation  of 
Caputo  et  al.  (2005,  their  Table  2).  We  again  fundamentalized 
the  pulsation  period  of  Polaris,  and  using  its  intensity-averaged 
value  of  Mv  we  find  Mp( PLC)  =  5.1  ±  0.4  M©.  The  pulsation 
mass  of  Polaris  can  also  be  estimated  using  the  predicted  period- 
mass-radius  (PMR)  relation  for  FO  Cepheids  of  Bono  et  al. 
(2001a),  along  with  the  radius  of  Polaris,  R  =  46  d=  3  RQ, 
measured  interferometrically  by  Nordgren  et  al.  (2000).  This 
gives  Mp (PMR)  =  4.9  ±  0.4  M0. 

The  lower  pulsation  masses,  taken  at  face  value,  are  thus  in 
better  agreement  with  the  nominal  dynamical  mass  than  are  the 
higher  evolutionary  masses.  However,  the  current  la  range  of 
the  measured  dynamical  mass,  3. 1-6.7  MQ,  encompasses  the 
entire  range  of  theoretical  masses.  Thus,  our  discussion  serves 
mainly  to  emphasize  the  crucial  importance  of  reducing  the 
error  bars  through  continued  HST  high-resolution  imaging  of  the 
Polaris  Aa,  Ab  system.  In  addition,  our  companion  HST  FGS 
astrometric  program  will  provide  an  improved  trigonometric 
parallax.  No  doubt  that  there  will  also  be  future  improvements 
in  the  spectroscopic  orbit  (e.g.,  Turner  et  al.  2006;  Bruntt  et  al. 
2008).  Simulations  suggest  that  we  can  reduce  the  uncertainty  on 
the  Cepheid’s  dynamical  mass  to  below  ±0.6  M0.  This  would 
provide  a  major  constraint  on  the  evolution  of  intermediate-mass 
stars  and  the  physics  of  Cepheid  pulsation. 

6.3.  Issues  in  the  Evolution  of  Intermediate-Mass  Stars 

As  a  further  elaboration  of  the  importance  of  accurate 
dynamical  masses  for  Cepheid  variables,  we  summarize  the 
major  open  questions  in  the  calculation  of  evolutionary  tracks 
of  intermediate-mass  stars.  A  more  complete  discussion  is 
provided  by  Bono  et  al.  (2006). 

The  luminosity  of  an  evolved  intermediate-mass  star  is  related 
to  the  mass  of  the  He-buming  core.  The  physical  mechanisms 
affecting  the  core  mass  include  the  following. 

1.  “Extra-mixing”  of  hydrogen  into  the  core  through  convec¬ 
tive  core  overshooting  during  the  central  hydrogen-burning 
phases. 

2.  Mass  loss,  leading  to  a  lower  total  stellar  mass  at  the  same 
luminosity. 

3 .  Rotation:  the  shear  layer  located  at  the  interface  between  the 
convective  and  radiative  regions  enhances  internal  mixing, 
producing  a  larger  He  core  mass. 

4.  Radiative  opacity:  an  increase  in  stellar  opacity  causes  an 
increase  in  the  central  temperature,  enhanced  efficiency  of 
central  H-burning,  and  a  higher  core  mass. 

Here  we  briefly  discuss  a  few  recent  results  that  bear  on 
convective  overshoot  and  mass  loss. 

The  discussion  of  the  Polaris  mass  in  the  previous  subsection 
showed  that  the  inclusion  of  noncanonical  overshoot  gave  a 
better  agreement  with  our  preliminary  dynamical  mass.  This 
is  borne  out  by  a  mass  measurement  for  the  longer-period 
Cepheid  S  Muscae,  based  on  HST  Goddard  High  Resolution 
Spectrograph  (GHRS)  radial  velocities  of  its  hot  companion, 
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and  an  assumed  companion  mass  based  on  its  Far  Ultraviolet 
Spectroscopic  Explorer  (FUSE)  spectrum  (Evans  et  al.  2006). 
The  implied  mass  of  S  Mus  clearly  favors  mild  convective 
overshoot. 

For  mass  loss,  we  note  that  the  evolutionary  calculations 
discussed  in  the  previous  subsection  included  semiempirical 
mass-loss  rates  (Reimers  1975;  Nieuwenhuijzen  &  de  Jager 
1990).  These  rates  are  insufficient  to  resolve  the  discrepancy  be¬ 
tween  evolutionary  and  pulsational  masses.  However,  a  variety 
of  mostly  recent  observational  information  suggests  that  mass 
loss  from  Cepheids  may  be  significant.  At  least  two  Cepheids, 
SU  Cas  and  RS  Pup,  are  associated  with  optical  reflection 
nebulae  (see  Kervella  et  al.  2008  and  references  therein)  that 
may  represent  mass  ejection  from  the  Cepheids.  Moreover,  a 
large  circumstellar  envelope  around  the  Cepheid  i  Car  has  been 
detected  recently  by  Kervella  et  al.  (2006),  using  mid-infrared 
data  collected  with  the  mid-infrared  instrument  (MIDI)  on  the 
Very  Large  Telescope  Interferometer  (VLTI).  Circumstellar  ma¬ 
terial  has  also  been  detected  around  8  Cephei  and  Polaris  itself 
by  Merand  et  al.  (2006).  Recent  Spitzer  observations  of  a  sam¬ 
ple  of  Cepheids  (Evans  et  al.  2007)  have  likewise  revealed  an 
infrared  excess  in  the  direction  of  8  Cephei. 

On  the  theoretical  side,  a  recent  investigation  (Neilson  & 
Lester  2008)  indicated  that  the  coupling  between  radiative  line 
driving  (Castor  et  al.  1975)  and  the  momentum  input  of  both 
radial  pulsation  and  shocks  can  provide  mass-loss  rates  for 
Galactic  Cepheids  ranging  from  10“ 10  to  10-7  M©  yr-1.  This 
finding,  together  with  typical  evolutionary  lifetimes  (e.g.,  Bono 
et  al.  2000,  Table  7),  indicates  that  classical  Cepheids  may  in 
fact  be  capable  of  losing  10-20%  of  their  mass  that  would  be 
needed  to  resolve  the  discrepancy. 

7.  SUMMARY 

The  results  of  this  study  are  as  follows. 

1.  We  have  used  UV  imaging  with  the  ACS/HRC  onboard 
the  HST  to  make  the  first  direct  detection  of  the  close 
companion  of  the  classical  Cepheid  Polaris. 

2.  We  confirm  orbital  motion  in  a  retrograde  sense,  based  on 
two  observations  a  year  apart. 

3.  By  combining  our  HST  measurements  with  the  single- 
lined  spectroscopic  orbit  (Kamper  1996)  and  the  FK5  and 
Hipparcos  proper  motions  (Wielen  et  al.  2000),  we  derive  a 
dynamical  mass  for  the  Cepheid  Polaris  Aa  of  4.5^' 4  M© — 
the  first  purely  dynamical  mass  for  any  Cepheid. 

4.  The  dynamical  mass  is  smaller  than  values  estimated  from 
either  pulsational  properties  or  evolutionary  tracks,  but  the 
error  bars  are  still  large  enough  that  the  discrepancies  have 
not  achieved  statistical  significance. 

5.  The  close  companion  Polaris  Ab  has  a  dynamical  mass 
of  1.26^0  07  M©.  This  is  consistent  with  a  spectral  type  of 
about  F6  V,  inferred  from  the  UV  brightness  of  Ab. 

6.  The  more  distant  and  well-known  companion  Polaris  B  has 
a  UV  flux  consistent  with  its  known  spectral  type  of  F3  V, 
lying  at  the  same  distance  as  the  Cepheid.  The  proper 
motion  of  Polaris  B  is  shown  to  be  very  similar  to  that 
of  Aa,  Ab,  consistent  with  motion  in  a  wide  but  bound  orbit 
around  the  close  pair. 

7.  Continued  HST  imaging,  including  two  more  observations 
that  have  been  approved  for  our  own  program,  will  decrease 
the  errors  on  the  dynamical  mass  of  Polaris,  allowing  a 
critical  test  of  stellar-evolution  theory  and  the  influence  of 


such  effects  as  convective  overshoot,  mass  loss,  rotation, 
and  opacities. 
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